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correlation function tomography 
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Abstract. During the past few years, secure detections of cosmic shear have been obtained, manifest in the corre- 
lation of the observed ellipticities of galaxies. Constraints have already been placed on cosmological parameters, 
such as the normalisation of the matter power spectrum ag. One possible systematic contaminant of the lensing 
correlation signal arises from intrinsic galaxy alignment, which is still poorly constrained. Unlike lensing, intrinsic 
correlations only pertain to galaxies with small physical separations, the correlation length being a few Mpc. 
We present a new method that harnesses this property, and isolates the lensing and intrinsic components of the 
galaxy ellipticity correlation function using measurements between different redshift slices. The observed signal is 
approximated by a set of template functions, making no strong assumptions about the amplitude or correlation 
length of any intrinsic alignment. We also show that the near-degeneracy between the matter density parameter 
f2 m and as can be lifted using correlation function tomography, even in the presence of an intrinsic alignment 
signal. 
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The tidal gravitational field of mass inhomogeneitics dis- 
torts the images of distant galaxies, resulting in correla- 
tions in their observed ellipticities. This cosmological weak 
lensing signal, or cosmic shear, depends upon cosmologi- 
cal parameters and the matter power spectrum (Blandford 
et al. 1991; Miralda-Escude 1991; Kaiser 1992). In 2000, 
four teams announced the first detections of cosmic shear 
(Bacon et al. 2000; Kaiser et al. 2000; van Waerbeke et al. 
2000; Wittman et al. 2000; Maoli et al. 2001), and more re- 
cently measurements at arcminute scales have been made 
using the HST (Hammcrlc et al. 2002; Refregier et al. 
2002). Interesting constraints have already been placed on 
the matter power spectrum normalisation as, and cosmic 
shear is also particularly sensitive to the matter density 
parameter Q m , the power spectrum shape parameter T 
and the source redshift distribution (e.g. van Waerbeke 
et al. 2002a; Hoekstra et al. 2002). Future multi-colour 
surveys will cover hundreds of square degrees, and have 
the potential to place tight constraints on cosmological 
parameters particularly when combined with results from 
the CMB, SNIa and galaxy surveys (Mellier et al. 2002; 
van Waerbeke et al. 2002b). For example, van Waerbeke 
et al. (2002a) compared their constraints on f2 m and as 
from lensing with those of Lahav et al. (2002) from the 
CMB, noting their near orthogonality. 



Various statistical measures of the cosmic shear have 
been suggested; here we focus on the two-point shear cor- 
relation function £ + (hereafter denoted by £), which is con- 
venient since it is insensitive to gaps in the data field, un- 
like integrated measures such as the aperture mass statis- 
tic A4 ap (e.g. Schneider et al. 1998). 

A possible systematic contaminant of the lensing cor- 
relation function £ L is intrinsic alignment, which may 
arise during the galaxy formation process. This has been 
subject to numerical, analytic and observational studies 
[e.g. Croft & Metzler 2000; Heavens et al. 2000 (HRH); 
Crittenden et al. 2001; Catelan et al. 2001; Mackey et al. 
2002; Brown et al. 2002; Jing 2002; Hui & Zhang 2002], 
where amplitude estimates span a few orders of magnitude 
due to differences in the mechanism assumed to be re- 
sponsible, and the type of galaxy considered. Nevertheless, 
these studies agree that the intrinsic correlation signal 
can dominate the lensing signal for surveys with (z) 0.5. 

Any correlation in ellipticities due to intrinsic align- 
ment only arises from physically close galaxy pairs, 
whereas £ L is sensitive to the integrated effect of the den- 
sity fluctuations out to the redshift of the nearer galaxy. 
It has been shown that photometric redshift information 
could be used to suppress £}, by downweighting or ignor- 
ing galaxy pairs at approximately the same redshift (King 
& Schneider 2002), or by downweighting nearby pairs and 
subtracting a model of the intrinsic alignment signal from 
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the observed ellipticity correlation function (Heymans & 
Heavens 2002). 

Motivated by the fact that intrinsic galaxy alignment 
is not yet well understood, we present a new method to 
isolate the intrinsic and lensing-induced components of 
the galaxy ellipticity correlation function. This method 
assumes that photometric redshift information is avail- 
able, so that the correlation function can be measured be- 
tween different redshift slices. However, no specific model 
for intrinsic ellipticity correlation (for instance its corre- 
lation length or redshift evolution) needs to be adopted. 
In the next section we outline the method and in Sect. 3 
we present some results in the context of a possible future 
survey. We discuss the results in Sect. 4. 



2. General method and specific assumptions 

In this section we outline a method to separate and ex- 
tract the intrinsic and lensing components of the galaxy 
ellipticity correlation function. The general method is de- 
scribed in Sect. 2.1, and in Sect.|2.2| and Sect. 2.3 we state 



the assumptions particular to this work. 



2.1. Method 

The ellipticity correlation function for galaxies with angu- 
lar separation 6 and at true redshifts Zi, Zj, is composed 
of a lensing and an intrinsic signal 



£(0, Zi, Zj ) = £ L (0, z h Zj) + £ l (9, Zi, Zj ) 



(1) 



As noted in Sect. 1, the origin and behaviour of is 
not yet well understood. £ L (0, Zi,Zj) is related to the 3- 
dimcnsional matter power spectrum P$ through 
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where Hq and f2 m are the values of the Hubble parameter 
and matter density parameter at the present epoch, and 
a(w) is the scale factor at comoving distance w, normalised 
such that a(0) = 1 today. Jo is the 0-th order Bessel func- 
tion of the first kind, and £ is the angular wave-vector. 
We denote the comoving distance at by Wi, and the 
function f(w) is the comoving angular diameter distance, 
which depends on the spatial curvature K: 

f(w) = \w (K = 0) } . (3) 

[ (-A')- 1/2 sinh[(-X) 1 /2 w ] (K<0) J 

The function R(w,w') = f(w' — w)/f(w') is the ratio of 
the angular diameter distance of a source at comoving 
distance w' seen from a distance w, to that seen from 
w = 0. 



Next, we account for the availability of photometric 
redshift estimates rather than spectroscopic ones. The 
galaxy ellipticity correlation function becomes 

i;{0,Zi,Zj) = J dzi J dzjp(zi,Zj\zi,Zj,6)£(6,zuZj) , (4) 

where p(zi, Zj\z~i, Zj, 9) is the probability to have true red- 
shifts Zi and Zj given photometric estimates z, and Zj. 
This is given by 

p(zi,Zj\zi,z~j,6) = (5) 

p(zj\zj) p(zj\zj) [1 + g gg (r)] 
/ dzi J dzj p{z l \z i )p{z i \zj) [1 + fgg(r)] ' 

where £ gg is the galaxy spatial correlation function, which 
may also include redshift dependent evolution, and r is 
the comoving separation of the galaxies in a pair. 

We now assume that £(0, Zi,z~j) is available on a 3- 
dimensional grid of Nk angular separation bins of width 
A8 centred on 0k (index K), and Nz photometric redshift 
bins of width Az centred on each of z\ (index /) and Zj 
(index J). This could either correspond to an observed sig- 
nal £fj K , or to a theoretical prediction (£ mod ) IJK which 
we want to compare with the observed signal. 

We assume that both the lensing and intrinsic correla- 
tions can be written in terms of sets of template functions 
A,, and B„ 



£ h (9,Zi,Zj) = y^a n A n (6,Zi,Zj) , 



(6) 



,) = ^Tb n B n ( 



n=l 



where a n and b n are the amplitudes of the n-th lensing and 
n-th intrinsic template functions. The template functions 
are fairly arbitrary, and extra functions can be added as 
required, to span the range of pla usib le models. We de- 
scribe our choice of models in Sect. 2.3 below. 



A suitably chosen single index m identifies correlations 
between bins with redshift indices I, J and angular sep- 
aration index K. In total, there are Nm = Nz {Nz + 
1)Nk/2 such independent measurements. The total of 
N = Nl + Nj gridded template models for the correla- 
tion functions (Ai...Aj^ L , Si...Bjv 7 ) can be written as an 
Nm x N so-called design matrix M. , and their amplitudes 
(ai...aif L ,bi...bff t ) as an A-dimcnsional column vector Q 
so that 



<r od ) r 



(7) 



Our aim is to recover in terms of the template func- 
tions. Using the method of least squares, the best-fit esti- 
mates of the Q n are those values Q n which minimise 

S = (r bs - MG) C- 1 {£ ohs - MQ) , (8) 

where C is the covariance matrix. Since (£ mod ) m is a linear 
combination of template functions, the linear least squares 
estimators are 



g = (m t c~ 1 m) 1 m t c- 1 $, 



1 ^obs 



(9) 
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After obtaining C, combined with an observed correlation 
function £ obs and a design matrix M, we have the neces- 
sary machinery to obtain Q i.e. the projection of £ obs into 
the template functions. This means that separate fits for 
the lensing and intrinsic contributions can be obtained. 

2.2. Covariance matrix 

To evaluate (J$), we need the covariance matrix C mm t. 
The covariance matrix could be calculated using the 
method described in Schneider et al. (2002), where C is 
expressed as integrals over (products of) correlation func- 
tions. However, since the method presented here would re- 
quire the calculation of very many elements of the covari- 
ance matrix, owing to the redshift slicing, we decided to 
use, as a first step, a simplified model for C. This consists 
of neglecting the cosmic variance contribution to C, and 
thus consideration of the (diagonal) elements of C coming 
from the intrinsic ellipticity dispersion of the source galax- 
ies. This is the dominant contribution to the covariance 
at small angular scales; at larger angular scales, the cos- 
mic variance terms start to dominate, with the transition 
angular scale depending on the survey geometry (Kaiser 
1998; Schneider et al. 2002). Here, we consider inde- 
pendent fields, and take Mr = 300. Therefore, we expect 
the cosmic variance not to be very much larger than the 
intrinsic ellipticity noise on the angular scales considered. 

With this approximation, the elements of the covari- 
ance matrix are 

2 (<r 6 2 /2) 2 

Cmm< = 77 , , [5 m ,rn> (1 + fa I'Sj J>5l j) , (10) 

A/p(m) 

where there is an extra contribution from auto-variance 
terms. The elements of the inverse covariance matrix in 
this case are simply (C~ 1 ) mm , = S mm >/C mm . The galaxy 
ellipticity dispersion is denoted by a e . In bin m, the num- 
ber of pairs is given by 

Af p (m)=Af F [n p{zi)Az] [n p(zj)Az] L 4 —t(- \ ,(11) 

where n is the galaxy number density, p(zj) is the redshift 
probability density for redshift bin /, and L is the extent 
of the field, assumed to be square. t(9/L) is a function 
that takes into account the fact that fields have finite ex- 
tent, and it must be evaluated numerically. Fig.^] shows 
the function t(8/L); note the limiting case where pairs 
separated by more than \/2L do not occur. 

2.3. Basis models 

In order to illustrate the general method described above, 
we will choose a simple, restricted set of template func- 
tions which share the approximate functional behaviour 
expected from the real correlation functions, both intrin- 
sic and lensing. For the latter, we simply take a small 
number of CDM cosmologies and consider their correla- 
tion functions as a template set. For the former, simple 
exponentials with redshift dependence are chosen. 




0.2 0.4 0.6 0.8 1 1.2 1.4 

9/L 

Fig. 1. The function r(6/L) (as defined in the text) which 
accounts for the finite extent of the data field when calcu- 
lating the number of pairs in a given bin. 

The lensing template functions A n used here are the 
gridded ^(0, Zi, Zj) for 3 models of the underlying cosmol- 
ogy: ACDM (fl m = 0.3, Q A = 0.7), OCDM (fi ro = 0.3, 
il A = 0) and crCDM (tt ni = I, Q A = 0); for simplicity, 
(j$ = 0.9 and T = 0.21 for each model. We use the Bardeen 
et al. (1986) transfer function to describe the evolution of 
the 3-dimensional power spectrum, along with the pre- 
scription of Peacock & Dodds (1996) for evolution in the 
non-linear regime. The required lensing correlation func- 
tions are calculated using the relationship between the 
power spectrum and £ L given in (ph, and then integrated 
over the photometric redshift uncertainties as in (||) . Here 
it is assumed that p(z\z) is a Gaussian with dispersion 
Cphot, centred on z. 

Nine template models B n for the intrinsic alignments 
are considered. First, the true spatial intrinsic correlation 
function is parameterised in terms of a correlation length 
i?corr and an exponent a: 

T?(r, z) = {\ + z av ) Q [exp (-r/iZoorr)] , (12) 

where z av is the mean redshift of galaxies in a pair and 
r is their comoving separation. We use the approxima- 
tion r 2 = (w l — Wj) 2 + 9 2 f 2 [(wi + Wj)/2\. R COYI was 
taken to be [1, 3, 10]ft. _1 Mpc and a to be [— 1, 0, 1]. 
The availability of photometric redshift estimates is then 
accounted for by integrating rj(r, z) as in (Q), and finally 
we obtain each of the model correlation functions on a 
grid. Note that the intrinsic models are calculated us- 
ing the relationships for the distances f(w) pertaining to 
the ACDM cosmology described above, and that we use 
CggW - (r/bh^Mpc)- 1 - 8 . 

3. Results 

The results are presented in the context of a possible fu- 
ture multi-colour cosmic shear survey, with a field size 
L = 14', and A/jr = 300 independent pointings (i.e. the 
largest scale on which the ellipticity correlation function 
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Fig. 2. The lensing and intrinsic cross-correlation functions measured between three pairs of photometric redshift bins 
for the ACDM cosmology and HRH intrinsic alignment model. Open (filled) squares correspond to the true lensing 
(intrinsic) signal for bins centred on z% — 1.07 and Zi = 1.13. Open (filled) triangles show the true lensing (intrinsic) 
signal for bins centred on Zi — 0.62 and Zi = 0.68. Open (filled) inverted triangles show the true lensing (intrinsic) 
signal between bins centred on Zi — 0.62 and Z{ — 1.13; since the intrinsic signal is so low, this is plotted on the inlay 
panel along with the lensing signal. The sets of lines associated with each of the true model symbols are the recovered 
best-fit lensing (intrinsic) cross-correlation functions to noisy realisations of for the pairs of photometric redshift 
bins, using the procedure described in the text. 
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Fig. 3. The lensing correlation functions for the three cosmologies used for the construction of template functions, 
and for the other cosmologies considered in Sect. 3.1, plotted for the bins centred on Zi — 1.07 and Zi — 1.13. Filled 
squares correspond to crCDM, the filled pentagons indicate ACDM and filled diamonds correspond to OCDM. The 
dashed and solid lines indicate the flat cosmologies with (i) cr 8 = 0.71, fl m = 0.33 and T = 0.215, and (ii) og = 0.8, 
O m = 0.5 and T = 0.3, respectively. 



Lindsay King & Peter Schneider: Galaxy ellipticity correlation function tomography 



5 




5.0E-04 



0.0E+00 



6 (arcmin) 

Fig. 4. The lensing and intrinsic correlation functions between different redshift bins for the ACDM cosmology and 
Jing (2002) intrinsic alignment model. The symbols and lines have the same meaning as in Fig.0. 
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Fig. 5. The lensing and intrinsic correlation functions between different redshift bins for a flat cosmology with O m = 
0.33, ag = 0.71 and T = 0.215 and an intrinsic alignment model from HRH. The symbols and lines have the same 
meaning as in Fig.||. 



is available is \[2 x 14'). A galaxy number density of 
30 arcmin~ 2 and ellipticity dispersion of er e = 0.3 are used 
throughout. The value of er p hot = 0.1 is chosen, typical of 
that obtained with current SED fitting procedures such as 
hyper-Z using a wide range of optical and near-infrared fil- 
ters (Bolzonella et al. 2000). There are N z = 65 redshift 
slices between z = 0.2 and 2.12, and Nr = 25 angular 
separation bins between 0'.3 and 15'. 



The galaxy redshift distribution follows the parame- 
terisation suggested by Smail et al. (1995), i.e. p(z) = 
(3/[z o T (3/(3)](z/za) 2 exp(-(z/zo)P), where 1^ denotes 
the gamma function. We take j3 = 3/2 and z = 2/3 
yielding (z) « 1. 

Two applications of the technique are presented. We 
start by asking how well the intrinsic and lensing contri- 
butions to the galaxy ellipticity correlation function can 
be separated, using the information contained in correla- 
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Fig. 6. The lensing and intrinsic correlation functions between different redshift bins for a flat cosmology with f2 m = 
0.5, eg = 0.8 and T — 0.3 and an intrinsic alignment model from HRH. The symbols and lines have the same meaning 
as in Fig. |[ 



tion functions between different redshift slices. Secondly, 
it has been shown that the degeneracy between f2 m and 
us can be partly lifted when redshift estimates for source 
galaxies are available (e.g. van Waerbeke et al. 2002a). 
We illustrate the use of correlation function tomography 
in this respect. 



3.1. Isolating the intrinsic correlation signal 

Now we consider input "observed" correlation functions 
comprising a lensing and an intrinsic contribution, to 
see how well the individual signals are recovered in terms 
of template functions using (^) . First it was checked that 
in the absence of noise, Q n = Q n , when £ obs is composed 
of a lensing and an intrinsic model contained in the set of 
templates. 

The intrinsic alignment model for spirals from HRH, 
r)(r, z) — 0.012 exp(— r/1.5 ft. _1 Mpc) was then used to 
obtain £ J , and £ L was calculated for a ACDM cosmol- 
ogy. Random gaussian distributed errors with dispersion 
a = C^;^j were added to these correlation functions giving 
noise realisations, and best-fit parameters Q n were recov- 
ered for each of these. Fig. |^ shows the (noise-free) input 
and recovered intrinsic and lensing correlation functions 
between three combinations of redshift slices: one close 
pair at z ~ 1, one close pair at z ~ 0.6, and the corre- 
lations between slices at z ~ 1 and z ~ 0.6. The lens- 
ing correlation function for each of the three cosmologies 
used in the construction of template functions is shown 
in Fig.^, plotted for the slices at z ~ 1.0. The intrin- 
sic correlation signal surpasses the lensing signal out to 
several arcminutes for both the low-redshift and the high- 



redshift bins. Considering bins with a large separation in 
redshift reduces the intrinsic signal to a negligible level, as 
expected. Even with our limited set of template functions, 
the reduced x 2 values of the recovered fits to the noise re- 
alisations are w 1. Also note that the intrinsic signal can 
be well represented in terms of the template functions, 
although it is not contained in the template set. 

Since the current template set for the intrinsic align- 
ment signal contains exponentially-decaying models with 
different scale-lengths, we now consider how well the 
method fares if the true signal is a power-law instead. 
The intrinsic alignment model is taken from Jing (2002); 
we use rj(r, z) = 0.288/[r°' 4 (7.5 1 ' 7 +r 1 - 7 )], with r measured 
in units of h~ l Mpc. Again, the true cosmology is ACDM. 
Noise was added and best-fit parameters recovered in the 
same manner as described above. Fig. || shows the (noise- 
free) input and recovered intrinsic and lensing correlation 
functions between the same three combinations of redshift 
slices as for Fig.||. Even though the functional form of the 
true intrinsic signal is quite different from the template 
models, the best-fit intrinsic models are still rather close 
to the noise-free model and more importantly, the lensing 
signal is again well recovered. To assess the difference be- 
tween using the Jing model rather than the HRH model 
for intrinsic alignments, the reduced x 2 values of the re- 
covered fits were determined for 1000 noise realisations of 
each, keeping the same ACDM cosmology. The reduced \ 2 
value is lower for the HRH realisation in nearly all cases, 
since this is an exponentially-decaying model for which 
the templates are better adapted. The x 2 values for the 
HRH realisations closely follow the theoretically expected 
distribution. Although the intrinsic models can be distin- 
guished statistically, the difference in the mean x 2 values 
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of the two sets of 1000 realisations is only « 15% of their 
dispersion. In practice, several families of functional forms 
could be taken for the template set. 

The next example for £ obs again uses the HRH model 
for intrinsic alignments, but this time a different flat cos- 
mology with us = 0.71, f2 m = 0.33, T = 0.215 was used 
for the lensing correlations; the lensing signal for this cos- 
mology was not part of the template set. Fig.|| shows the 
lensing correlation function for this cosmology, plotted for 
the the slices at z « 1.0. Noise was added using the same 
random seeds as above, and best-fit parameters recovered 
as before. Fig.|5| shows results for the same combination 
of redshift slices as for Fig.^J. The power spectrum corre- 
sponding to this cosmology has a lower normalisation - 
the ~ 20% reduction in cr 8 is not offset by the 10% in- 
crease in f2 m ; this is evident in the lower amplitude of the 
lensing correlation functions. Again, the recovered param- 
eters are consistent with the true "observed" correlation 
functions and the reduced y 2 values ss 1. 

Whereas the foregoing cosmological model was quite 
different in amplitude of the power spectrum we now con- 
sider a model which also differs in the shape of the power 
spectrum, to test the robustness of our method. Keeping 
the HRH model for intrinsic alignments, a flat cosmol- 
ogy with erg = 0.8, ft m = 0.5, r = 0.3 was next used 
for the lensing correlations. Fig.|^ shows the lensing cor- 
relation function for this cosmology, plotted for the the 
slices at z ~ 1.0. Noise was added and best-fit param- 
eters recovered as described above. Fig.|| shows results 
for the same combination of redshift slices as for Fig.g. 
The lensing signal is again well represented in terms of 
the basis functions, even though it is rather different to 
any of the cosmology templates. Hence, despite the fact 
that our set of template functions is quite restrictive, we 
have demonstrated that it provides enough flexibility to 
provide accurate fits to the correlation functions of quite 
different comological models. 

3.2. Breaking the il m -ag degeneracy 

The near degeneracy, in the absence of redshift informa- 
tion, between f2 m and as for two example flat cosmolo- 
gies is illustrated in Fig.fj], taken from King & Schneider 
(2002; where details for its calculation can be found). 
One cosmology is the fiducial ACDM model: Q m = 0.3, 
erg = 0.9, and the other is an almost degenerate model 
with Q m = 0.4 and as — 0.78. To obtain these correlation 
functions, we assume the same prescription for the power 



spectrum outlined in Sect. 2.3 above. The source popula- 
tion has a redshift probability distribution with (z) = 1, 
but with no individual photometric redshift estimates as- 
sumed to be available. 

In order to see how well the fiducial and degenerate 
models could be distinguished using correlation function 
tomography, we took the gridded correlation functions for 
each model and added each of these to the gridded HRH 
intrinsic correlation function, giving 

(£ m ° d )m for each - For 




(arcmin) 



Fig. 7. This figure shows the lensing correlation functions 
for our fiducial ACDM model: tt m = 0.3, as — 0.9 (dashed 
line), and for a degenerate model with fi m = 0.4, as = 
0.78 (solid line). Both functions were predi cted using the 
matter power spectrum described in Sect. 2J3, and for a 
redshift distribution with (z) = 1. 



simplicity, we refer to these combinations as £acdm an d 
£dcgcn- A first set of simulations involved using a set of 
ten template functions containing the nine models for in- 
trinsic alignments, along with the ACDM lensing model 
for the lensing template. In turn, 1000 noise realisations 
of £acdm an d £dcgcn wer e generated using the same ran- 
dom seeds in both cases, and the best-fit amplitudes Q n 
for the template functions recovered. The process was 
repeated, this time using the gridded degenerate model 
as the lensing template function, in place of the fiducial 
ACDM model. The histograms of (i) Ax 2 (ACDM-degen) 
and (ii) A% 2 (degen — ACDM), corresponding to the dif- 
ference in goodness-of-fit for the noise realisations when 
€ acdm an( l £degen are used in the template set, are shown 
in Fig.||. When the fiducial ACDM model (degenerate 
model) is the best-fit and is contained in the template 
set, values of A% 2 in Fig.|| should be negative. This gives 
a measure of our ability to differentiate between models 
using correlation functions between redshift slices. In the 
first histogram, when the model for £acdm IS contained 
in the template set, in 95.6% of cases the noisy ACDM 
correlation functions are better fit. Also, when £<j egen is in 
the template set, 96% of the noisy degenerate correlation 
functions have better fits. Hence, within the assumptions 
we made and in the presence of an intrinsic alignment sig- 
nal, these two cosmological models could be distinguished 
at the ~ 2cr-level. 

The foregoing example for distinguishing between two 
cosmological models should serve as a general illustration 
only. Owing to our simplified ansatz for the covariance ma- 
trix C, which contains intrinsic ellipticity dispersion only, 
a more detailed investigation is not warranted here. In 
practice, cosmic variance would need to be taken into ac- 
count in the covariance matrix when realistic constraints 
on cosmological parameters are to be derived. We are cur- 
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Fig. 8. These histograms show the frequency (NR) of val- 
ues of Ax 2 for two sets of 1000 noise realisations. The first 
(second) set of realisations indicated by the solid (dashed) 
histogram has the fiducial (degenerate) ACDM cosmol- 
ogy in the template set and A% 2 is calculated between 
the best-fit ACDM (degenerate) as opposed to degenerate 
(ACDM) recovered model. 

rently investigating ways to obtain a far more realistic rep- 
resentation of the covariance matrix, to be used in a study 
of the accuracy of cosmological parameter determination. 

4. Discussion and conclusions 

It has been suggested that the lensing correlation func- 
tion may be contaminated by intrinsic galaxy alignments. 
Since cosmic shear probes the matter power spectrum and 
enables constraints to be placed on cosmological parame- 
ters such as cr§ and fi m (e.g. van Waerbeke et al. 2002a), 
it is vital to have the ability to isolate the contribution 
from intrinsic galaxy alignments in order to remove this 
systematic. Of course, intrinsic alignment is interesting in 
its own right: its amplitude as a function of physical sepa- 
ration and its evolution with redshift provides clues about 
the galaxy formation process. 

We have demonstrated that measuring galaxy elliptic- 
ity correlation functions between redshift slices would en- 
able the intrinsic and lensing contributions to be disentan- 
gled. The total signal is decomposed into template func- 
tions, and the fact that intrinsic alignments operate over 
a limited physical separation enables the intrinsic compo- 
nent to be isolated and subracted from the total signal. 
Our knowledge of the amplitude of intrinsic alignments 
is limited, but no strong assumption about the behaviour 
of the intrinsic alignment signal needs to be made. Here 
we considered a modest number of template functions, 
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Fig. 9. The ratio of partial derivatives of £ to £ for our 
fiducial ACDM model, with respect to parameters (i) f2 m , 
(ii) as, (iii) T and (iv) n. 



which can easily be augmented to cover a wider range 
of functional forms. For example, any intrinsic alignment 
signal arising at the epoch of galaxy formation may be 
suppressed by subsequent dynamical interaction, perhaps 
most pertinent to galaxy pairs with extremely small phys- 
ical separations. In fact, if the reduced x 2 of the best fit is 
significantly larger than 1, this indicates that additional 
template functions need to be included. 

Our choice of template functions is of course fairly 
arbitrary. We have taken functions which have approxi- 
mately the behaviour expected from a cosmic shear mea- 
surement. Alternatively, one could consider a set of generic 
basis functions, which however, owing to the dependence 
on three variables, would require a fairly large set of func- 
tions. Another natural choice of the template functions 
could be the following: assuming a reasonable guess for 
the cosmological model, characterised by the parameters 
7To, the correlation function for neighbouring models could 
be written as 



(13) 



UTV 



and therefore, the set consisting of £(7To) and its partial 
derivatives with respect to the relevant cosmological pa- 
rameters would provide a useful set of template functions. 
As an illustration, in Fig.|| we have plotted the partial 
derivatives of £(#) for the case where no redshift infor- 
mation is available, i.e. the derivatives of the redshift- 
averaged correlation function, again with (z) ~ 1. For 
our cosmological model (characterised by 7Tq) we take the 
fiducial ACDM model. Derivatives are taken with respect 
to (i) f2 m , (ii) as, (iii) T and (iv) the primordial spec- 
tral index n (where our fiducial model is scale-invariant 
n = 1). We plot the ratio where the numerator is 

denoted by A(£, 7rj). In the limiting case where one curve 
is a scaled version of another, A(£,7r,) cx A(£,7Tj), it is 
impossible to first order to break the degeneracy between 
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parameters 7Tj and ttj. We again see a nice illustration of 
the near degeneracy between f2 m and a%, manifest in the 
similarity between A(£, Sl m ) and A(£, cr 8 ). Since erf enters 
into the linear power spectrum just as a pref actor, on large 
angular scales (i.e. in the linear regime), the curve for as 
tends to a constant, A(£,er 8 )/£ — * 2(7^ 1 for large sepa- 
rations. Another feature to note is that on the scale of a 
few arcminutes the curves for T and n change sign, imply- 
ing that there is less degeneracy between these parameters 
and either of f2 m or as- 

There are several ways in which the method and re- 
sults discussed here could be used. One way would be to 
consider the resulting split into intrinsic correlations and 
lensing signal as the final result, and to compare the re- 
sulting functions with theories of galaxies formation which 
predict the intrinsic alignment signal, and cosmological 
models predicting the shear correlation function. The re- 
sulting fits are, however, difficult to interpret statistically, 
i.e. the error bars on the shear correlation function are 
difficult to obtain. An alternative would be to consider 
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the fitted intrinsic signal only, subtract it from the cllip- 
ticity correlation function, and consider the result as the 
shear correlation function, together with the correspond- 
ing error bars. Subsequently, the correlation function can 
then be used for the redshift-weighting method of King & 
Schneider (2002), of course yielding much smaller contri- 
butions from the intrinsic correlations than for the unsub- 
tracted data. Furthermore, the resulting model for the in- 
trinsic correlation function could also be used as input for 
the subtraction method discussed in Heymans & Heavens 
(2002). 

In addition to providing a key to the suppression of any 
intrinsic alignment signal, photometric redshift estimates 
enable much tighter constraints to be placed on cosmo- 
logical parameters obtained from cosmic shear surveys, as 
demonstrated by Hu (1999). Although our prime goal in 
this paper is not the constraint of cosmological parame- 
ters, we have illustrated that the degeneracy between fi m 
and as can be lifted by observing correlation functions 
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betweeiti redshift slices, even when an intrinsic alignment 
systematic is present. 
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